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I calculate the kaon B-parameter Bk, defined via §(t)1k fx) 2 Bk = {K\syp,(l — 7s)ds7 M (l — 
'j5)d\K}, with a lattice simulation in quenched approximation. The lattice simulation uses an action 
possessing exact lattice chiral symmetry, an overlap action. Computations are performed at two 
lattice spacings, about 0.13 and 0.09 fm (parameterized by Wilson gauge action couplings /3 = 5.9 
and 6.1) with nearly the same physical volumes and quark masses. I describe particular potential 
difficulties which arise due to the use of such a lattice action in finite volume. My results are 
consistent with other recent lattice determinations using domain- wall fermions. 

o ■ 

O ' !• INTRODUCTION 

<n : 

The kaon B-parameter Bk, defined as ^{rriKfK) 2 BK = (K\sjfj,(l — 75)^57^(1 — j^dlK), is an important ingredient 
in the testing of the unitarity of the Cabibbo-Kobayashi-Maskawa matrix p| . It has been a target of lattice calculations 
since the earliest days of numerical simulations of QCD. Lattice calculations of Bk require actions with good chiral 
properties, since the matrix element of the four-fermion operator scales like the square of the pseudoscalar meson 
mass as that mass vanishes. If the lattice action does not respect chiral symmetry, the desired operator will mix with 
operators of opposite chirality. The matrix elements of these operators do not vanish at vanishing quark mass, and 
"JT therefore overwhelm the signal. 

There has been a continuous cycle of lattice calculations using fermions with ever better chiral properties. This 
calculation is yet another incremental upgrade, to the use of a lattice action with exact SU(Nf) ® SU(Nf) chiral 
symmetry, an overlap action. These actions have operator mixing identical to that of continuum-regulated QCD 

o : a 

The first lattice calculations of Bk were done with Wilson-type actions. Techniques for handling operator mixing 
have improved over the years, (for recent results, see 0-0) but this approach remains (in the author's opinion) 
arduous. 

Staggered fermions haw enough chiral symmetry at nonzero lattice spacing, that operator mixing is not 

a problem. One can obtain extremely precise values for lattice-regulated Bk at any fixed lattice spacing. However, 
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' to date, all calculations of Bk done with staggered fermions use "unimproved" (thin link, nearest-neighbor-only 
interactions), and scaling violations are seen to be large. For example, the JLQCD collaboration saw a thirty per 
cent variation in Bk over their range of lattice spacings. 
. ! Domain wall fermions pin chiral fermions to a four dimensional brane in a five dimensional world; chiral symmetry 
is exact as the length of the fifth dimension becomes infinite. For real-world simulations, the fifth dimension is finite 
and chiral symmetry remains approximate, though much improved in practice compared to Wilson-type fermions. 
Two groups fiol lllj have presented results for Bk with domain wall fermions. Ref. ^(J has data at two lattice 
spacings and sees only small scaling violations. There is a few standard deviation disagreement between the published 
results of the two groups. 

Finally, overlap actions have exact chiral symmetry at finite lattice spacing. All operator mixing is exactly as in 
the continuum Q . Two groups |l'2t ll.'3| have recently presented results for Bk using overlap actions, but the actions 
and techniques are completely different. The second of these is a preliminary version of the work described here. 

The lattice matrix elements must be converted to their continuum-regularized values and run to some fiducial scale. 
Matching coefficients can be computed perturbatively or nonperturbatively. For most standard discretizations of 
fermions, the matching factors ( "Z-factors" ) are quite different from unity. For these actions, perturbation theory 
is regarded as untrustworthy, and other methods must be employed. I, however, use an action in which the gauge 
connections are an average of a set of short range paths, specifically HYP-blocked links0|- ^ e f- US has computed 
the matching factors for operators relevant for this study, as well as the scale for evaluating the running coupling, using 
the Lepage-Mackenzie- Hornbostel 0,^3 criterion. The Z-factors are quite close to unity. A number of calculations of 
Z-factors for related actions 01 reveal that this behavior is generic for actions with similar kinds of gauge connections. 
In this work I compare perturbative and nonperturbative calculations of two Z-factors, for the local axial current and 
for the lattice-to- MS quark mass matching factor, and find reasonable agreement between them. 

Having exact chiral symmetry forces one to confront the problem that this calculation is performed in quenched 
approximation. These results will not be directly applicable to the real world of QCD with dynamical fermions. I 
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encounter difficulties in two places. First, there is no reason for the spectrum of quenched QCD to be identical to 
that of full QCD. Using different physical observables to extract a strange quark mass can (and does) lead to different 
values of this parameter. At small quark masses, my observed B parameter varies strongly with quark mass, and my 
prediction for Bk is sensitive to my choice of m s . 

The second difficulty is more fundamental. In order to extrapolate results to the chiral limit, one must use chiral 
perturbation theory. However, the symmetries of quenched QCD are different from the symmetries of full QCD. Not 
only are the leading coefficients expected to be different in quenched and full QCD, but the logarithmic contributions 
to any observable Q 



Q(mps) = A(l + fl-p- logm PS ) + 
J PS 



(1) 



can have different coefficients (different -B's in Eq. [Q, or different functional form (in the formula for m 
coefficient of log {m 2 PS ) is not mp s but a constant related to the quenched topological susceptibility 0,|2(J]). All of 
these differences are encountered in the analysis of the data. To produce a prediction for an experimental number 
using results of a quenched simulation involves uncontrolled, non-lattice-related phenomcnological assumptions. 

I conclude this introduction by presenting graphs which illustrate my results. Fig. ^shows results for Bk at various 
lattice spacings, for a selection of simulations which have reasonable statistics and small error bars. Fig. [21 presents 
results which are either extrapolated to the continuum limit, or presented by their authors as having small lattice 
spacing artifacts. 
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FIG. 1: Bk comparisons vs lattice spacing, from a selection of simulations with reasonably small error bars. Results are 
labeled diamonds and fancy diamond|!|, the fancy cross [T3| . octagons |T3| . the cross |Tll| . and squares (this work). 

In Section 2 I describe the action, simulation parameters, and data sets. Sec. 3 is devoted to a discussion of zero 
mode effects and my attempts to deal with them. Results relevant to Bk are presented in Sec. 4. In Sec. 5 I discuss 

3/2 3/2 

the chiral limit of Bk and of operators O r and O g , relevant for part of e'/e. My brief conclusions are given in 
Sec. 6. An Appendix compares perturbative and nonperturbative matching factors. 
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FIG. 2: Bk comparisons presented "as if" they were taken to the continuum limit. The label in parentheses characterizes 
the kind of lattice fermions used: W for Wilson, CI for Clover, KS for staggered, DW for domain wall, and OV for ove rlap 
fermions. References are GBS 0, JLQCD(W) 0, ContiS, KGS 0, JLQCD(KS) 0, CP-PACS |HJ, RBC [H, GGHLR (Tf, 
and MILC, this work. The points of Refs. Q, H, 0, US, 

and this work are the results of a a continuum extrapolation; all 

the rest are simulations at one lattice spacing. 



II. SIMULATION TECHNIQUES 
A. Data sets 



The data set used in this study is generated in the quenched approximation using the Wilson gauge action at 
couplings j3 = 5.9 (on a 12 3 x 36 site lattice), where I have an 80 lattice data set, and j3 — 6.1 (on a 16 3 x 48 site 
lattice) with 60 lattices. The nominal lattice spacings are a = 0.13 fm and 0.09 fm from the measured rho mass. 
Propagators for ten (/3 = 5.9) or nine {ft — 6.1) quark masses are constructed corresponding to pseudoscalar-to- vector 
meson mass ratios of mps/mv ranging from 0.4 to 0.85. The fermions have periodic boundary conditions in the 
spatial directions and anti-periodic temporal boundary conditions. I gauge fix to Coulomb gauge and take our sources 
to be Gaussians of size xo/a = 3, 4.125 at (3 = 5.9, 6.1 (where the quark source is $ = exp(— x 2 /x^)). 



B. Lattice action and simulation methodology 



The massless overlap Dirac operator is 



D(0)=x (l + 



(2) 
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where z = cI(—xo)/xq — (d—xo)/xo and d(m) = d + m is a massive "kernel" Dirac operator for mass m. The massive 
overlap Dirac operator is conventionally defined to be 

D(m q ) = (l-^-)D(Q) + m q (3) 

and it is also conventional to define the propagator so that the chiral modes at A = 2xq are projected out, 

D-\ mq ) = r. 1 -j-—(D-\m q ) - (4) 

l-m q /{2xo) 2x 

This also converts local currents into order a 2 improved operators |2l|. 

The overlap action used in these studies[22| is built from a kernel action with nearest and next-nearest neighbor cou- 
plings, and HYP-blocked links0|- HYP links fatten the gauge links without extending gauge-field-fermion couplings 
beyond a single hypercube. This improves the kernel's chiral properties without compromising locality. 

The "step function" (e(z) = z/Vz^z) is evaluated using the fourteenth-order Remes algorithm of Ref. Q (after 
removing the lowest 20 eigenmodes of z^z). This involves an inner (multimass j3|) conjugate gradient inversion step. 
It is convenient to monitor the norm of the step function |e(z)'i/'| 2 /|'i/'| 2 and adjust the conjugate gradient residue to 
produce a desired accuracy (typically 10~ 5 in e(z)ip). Doing so, I need about 16-18 inner conjugate gradient steps at 
(3 = 5.9 and 10-12 steps at (3 = 6.1. 

An important ingredient of this overlap program has been to precondition the quark propagator by projecting low 
eigenfunctions of the Dirac operator out of the source and including them exactly. This can in principle eliminate 
critical slowing down from the iterative calculation of the inversion of the Dirac operator. Of course, there is a cost: 
one must find the eigenmodes. My impression from the literature, plus my own experience, is that for the overlap with 
Wilson action kernel, this cost is prohibitive. However, my kernel action is designed to resemble the exact overlap 
well enough that its eigenvectors are good "seeds" for a calculation of eigenvectors of the exact action, and it is kept 
simple enough that finding its own eigenvectors is inexpensive. 

As a rough figure of merit, consider the 12 3 x 36, (3 — 5.9 data set. Computing the lowest twenty eigenmodes of 
the squared massless overlap Dirac operator takes about 8 time units, while the complete set of quark propagators 
from the lightest mass I studied to the heaviest takes about 16 time units times two (for two sets of propagators) per 
lattice. Fig. [3] shows the number of conjugate gradient steps needed to compute the quark propagator to some fixed 
accuracy, as a function of bare quark mass. If one would naively extrapolate the heavy quark results (for which low 
eigenmodes make only a small contribution) to lower masses, one would see a forty per cent reduction in the number 
of inverter steps needed at the smallest quark mass, vs an additional cost of 8/32 = 25 per cent per propagator set 
for the construction of eigenmodes. 



C. Correlation functions 



The "generic" four fermion operator one must consider is 

o = ) ® (#r 2 #>) (5) 

(the superscript labels flavor; the subscript, color). Special cases are (a) O = 0\. Y\ = T 2 = 7^(1 — 75), a = 5, (3 = 7; 
(b) O = O2: I"i = T 2 = 7 M (1 — 75), a — (3, 7 = 5; and (c) the isospin 3/2 operators for electroweak penguins, here 
written with the normalization conventions used in Refs. [llj and |25| 

°7 /2 = ^((^7p( 1 - 75)da)[(%37M(l + 75W) - (d/TY^O 1 +75)^)] + («o7At(l ~ l5)u a ) {u 0l ^ (1 + 75)^)) (6) 

and 

°8 /2 = \{{s a l^{l ~ 75)d/3)[(u/37 M (l +75M*) - (^7^(1 + 75)^)] + (ScO^l - 7 5 )W/3)(M/37 M (1 + 75K0)- (7) 

Bk is proportional to the matrix element of 0+ = 0\ + Oi- Because overlap fermions are chiral, one can extract 
Bk from the matrix element of the operator between zero momentum states. All the operators to be studied have 
only "figure-eight" topology matrix elements, where each field in the operator contracts against a field in the source 
or sink interpolating field. There are no penguin graphs (where fields in the operator contract against each other) as 
long as one works in the degenerate- mass limit. 
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FIG. 3: Number of conjugate gradient steps in the calculation of quark propagators, as a function of quark mass. The four 
curves are two different sources on each of two different lattices, and we are converging to a fractional accuracy of the squared 
residue of 1CF 14 for this test case. 

Broadly speaking, the matrix element of a four-fermion operator is computed by placing interpolating fields for a 
meson at two widely-separated locations on the lattice and contracting field variables between these sources and the 
operator, to construct an un-amputated correlator containing the operator. There are two commonly-used strategies 
for doing this: One possibility is to build all the quark propagators "at the operator" at one location on the lattice, 
and to join up pairs of propagators to make the mesons. The second method is to construct propagators from two 
well-separated sources and bring them together at the operator. The advantages of the first method are that one needs 
half as many propagators per lattice, and it is possible to project the whole calculation into a particular momentum 
eigenstate (by appropriately summing over locations of the meson interpolating fields). The major disadvantage of this 
method is that one is only measuring the matrix element of the operator at one location per propagator construction. 
With the second method one can average the location of the operator over all spatial and many temporal locations, a 
considerable gain in statistics. A disadvantage of the second method is that unless the source generates a hadronic state 
which is a momentum eigenstate the correlator will involve a mixture of the eigenstates. The dominant contribution 
to the correlator is from p = 0, but there will be a contamination from higher momentum states. I have chosen 
the second method, and will discuss below how I dealt with higher-momentum modes. I placed the two source time 
slices Nt/2 — 2 temporal sites apart; with the toroidal geometry there are two temporal regions where the operator is 
"between" the sources. I combine all two-point function data sets from the two sources in the fits. 

The correlator of two interpolating fields located at x, t = (0,0) and (0, T) with an operator S(x,t) summed over 
x is then 



c 3 (t,T) = j2(mT)s( X ,t)mo)) 



(8) 



Inserting complete sets of relativistically normalized momentum eigenstates, and assuming (h(pi)\S(p)\h(p2)) is pro- 
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portional to 5(p + p 2 — Pi), 

C 3 {t,T)^eM~E{mm\HM 2 ^^^{hmS\Kp)). (9) 

If the p = state dominates the sum, 

C 3 (t,T)=exp(-mT)K0|$|/i)| 2 -i- I i(/ i |5|ft). (10) 

(2m) z V 

However, "dominance in the sum" is controlled by the quantity exp((— E(p) — m)T). These sources do not make 
momentum eigenstates, and so thep = Bk signal is contaminated by ap ^ contribution. The minimum momentum 
in a finite box of side L is \p\ = 2ir/L. Because of the small masses involved, the contamination is not present at quark 
masses below 2-3 times the strange quark mass. It appears at bigger quark mass, because E{p) — m gets smaller as 
the pseudoscalar mass m grows. Fortunately, there are two inequivalent paths on the torus to disentangle the two 
"signals," and one can fit the Bk correlator to a sum of a p = term and a p — 2ir / ' L term. 



D. Fitting and Error Analysis 

In a Monte Carlo simulation different quantities measured on the same set of lattices are highly correlated, and it 
is preferable to analyze the data using covariant fits. However, in order to determine accurately the small eigenvalues 
of the correlation matrix, large statistical samples are needed. If the statistical sample is too small, these eigenvalues 
can be poorly determined, and the fits will become unstable. This can be a particular problem in fits which extract 
Bk, simply because they can involve many degrees of freedom (one three-point function and up to two two-point 
functions). My solution to this problem is to adopt the strategy of the analysis used in Ref. |2g. I compute the 
correlation matrix (the covariance matrix, but normalized by the standard deviation so that its diagonal entries are 
the identity), and find its eigenvalues and eigenvectors. I then reconstruct the matrix, discarding eigenvectors whose 
eigenvalues are smaller than some cutoff. This matrix is singular, so I reset the values of its diagonal entries to the 
identity again. I use this processed covariance matrix as an input to the fit. 

This is basically the strategy adopted in Ref. [2(j, except that here the cut on kept eigenvectors is chosen to be a 
fraction of the largest eigenvector, rather than a fixed value. I have taken the value of this cutoff to be 0.1, and have 
varied it to insure that fits are insensitive to it. 

(My experience in this work is that all fits to spectroscopic quantities, and fits to pairs of two point function, used, 
for example, to extract pseudoscalar decay constants, produce statistically identical results regardless of the number 
of eigenvalues discarded.) 

All chiral extrapolations of data are done by performing correlated fits to extract necessary physical parameters 
quark mass by quark mass, then performing extrapolations under a single-elimination jackknife. 

I have performed a number of different analyses of the data. Fits of the three point function with a particular 
sources (pseudoscalar or axial vector, labeled PS or A) can be combined with two-point functions with axial vector 
sinks to give the B-parameters directly, 

C 3 (f,T) = ^(Sexp(-mT) + C ie xp(-S(p)T)) (11) 

(E(p) 2 = m 2 + (2ir/L) 2 is the energy of the lowest nonzero momentum state) and 

C 2 (t, t') = Z(exp(-m(t - t') - exp(-m(N t - (t - t')). (12) 

This is done by performing a correlated fit with three (Z, B, mps) or four (add C\) parameters to C^{t, T), C2(t, 0), 
and C 2 (t,T). 

As a check, one can also simply perform a (jackknife-averaged) fit of the ratio C^{t, T)/[C 2 (t, 0)C 2 (t, T)} to a 
constant. The range of t is varied with a search for plateaus in the fit value. This method has been commonly used 
in many previous measurements of B parameters. When I compare this "ratio" fit to a full correlated fit to Bk, I 
always find consistency of the fitted B parameter. Generally, the error bar assigned to it by the ratio fit is about 
forty percent ot the size of the uncertainty of the correlated fit. However, I elect NOT to use these more optimistic 
estimations of the uncertainty because they do not attempt to include any of the correlations which are known to be 
present in the data. 
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One can also extract the matrix element of the operator by doing a correlated fit of C%(t, T) for a particular source 
and a two-point function which uses the same source as a sink. In this case 

C ^ T ) = T^ttA ) cxp(-mT) + C[ exp(-E(p)T)) (13) 
(2m) z V 

and 

Z' 

C 2 (t, t') = — (exp(-m(t - t') + exp(-m(T - (t - t')). (14) 
2m 

With these correlators, the B parameter can be extracted by doing a simultaneous fit to another correlator set, which 
gives fps and mps, and jackknife averaging. This increases the resulting uncertainty in Bk , so this method is not 
competitive with the three-propagator fits. 

At lower quark mass the size of the contamination falls to zero. Fits which include the C\ coefficient become 
unstable, because the fit is completely insensitive to it. The Hessian matrix develops a near-zero eigenvalue as a result 
of this insensitivity. Usually, this near zero mode just means that the C\ coefficient will be poorly determined, but a 
zero mode in a Hessian matrix is always dangerous, and can corrupt the whole fit. To determine the minimum quark 
mass where G\ is not needed, I performed fits where the inversion of the covariance matrix was done using singular 
value decomposition (clipping out eigenmodes with eigenmodes smaller than some minimum cut). I varied the ratio 
of smallest to largest eigenvalues kept in the 4x4 Hessian matrix. Any "reasonable" choice (conditioning number 
cut > 10~ 5 ) stabilizes the inversion by decoupling and freezing out any variation in C\. 

III. ZERO MODE EFFECTS 

Simulations with lattice fermions possessing an exact chiral symmetry done in finite volume have a "new" kind of 
finite volume artifact: the presence of exact zero modes of the Dirac operator, through which quarks can propagate. 
In some cases (the eta prime channel, for example) the zero modes contribute physics, but for the case of Bk they 
arc merely an annoyance which would disappear in the infinite volume limit. 

Things could be much worse. For example, for Wilson-type fermions the analogs of zero modes are configurations 
where the Dirac operator has real eigenvalues. Because these theories are not chiral, eigenvalues of the massless Dirac 
operator D are not "protected"- they can take on any real value. It can happen that this real value coincides with 
the negative of the simulation mass m. Then D + m is simply non-invertible. Typically, these eigenmodes occur at 
small quark mass, meaning that it is difficult (if not impossible) to do simulations there. 

Note that zero modes will also be present in simulations with nonzero mass dynamical fermions, although their 
effects arc likely to be less severe; the dynamical fermions suppress zero modes but (except at zero quark mass) do 
not eliminate them. 

There are basically five observables needed in a Bk calculation. 

• The pseudoscalar mass 

• The vector meson mass (to set the lattice spacing) 

• The matrix element of 0+ 

• The matrix element of the axial vector current 

• The pseudoscalar decay constant 

All can (in principle) be affected by finite volume zero modes artifacts at small quark mass. 

The most obvious way to check for finite volume effects is to do simulations at several volumes. I did not do that. 
One could use sources which do not couple to zero modes. This is the strategy pursued by Ref. |T^ : their meson 
interpolating fields are the linear combination 7075 — 70 • I did not do that either. The 70 part of the source couples 
only to the (heavy) scalar meson, not to the pseudoscalar meson. Its contribution to the pseudoscalar meson matrix 
element of the operator averages to zero but contributes noise in a finite data set. Because of the PCAC relation, the 
7075 source decouples for the pseudoscalar meson in the chiral limit, so that while one is decoupling zero modes, we 
are also decoupling from the desired channel. 

Instead, I used the fact that zero mode effects are different in different channels. I measured the same physical 
observable using different operators (presumably with different zero mode contributions) and looked for differences. 

When one fits correlators to extract pseudoscalar masses, one can perform the fits naively, without consideration 
of the zero modes. Then if zero modes are present in some channels, observables such as the pseudoscalar mass 
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inferred from these channels will show characteristic differences. As the zero mode correlator reflects some underlying 
structure of the gauge field responsible for configurations which support zero modes, it produces correlations on a 
scale which is independent of the quark mass. Thus, a plot of m 2 PS /m q would diverge at small m q inversely with m q 
(simply because mps would be independent of m q ). I searched for this suspicious behavior in plots of m 2 PS /m q vs 
m q (Fig. Q . Correlators include the pseudoscalar current (which has a contribution from the propagation of both 
quarks through zero modes), the pseudoscalar-scalar difference (in which all zero mode contributions cancel), and the 
axial current correlator (which has a "mixed" contribution with one quark propagating through zero modes and the 
other through nonzero modes: this contribution in independent of m q at small m q .) Within statistical uncertainties, 
I do not see any difference between these channels. I anticipate that this would not be the case at smaller quark mass 
or smaller volumes. Nevertheless, on theoretical grounds I elect to only extract pseudoscalar masses at small quark 
masses from the pseudoscalar-scalar difference. 
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FIG. 4: (a) (3 = 5.9 and (b) (3 = 6.1: (am PS ) 2 /{am q ) vs am q from pseudoscalar (crosses), axial vector (octagons) and 
pseudoscalar-scalar difference (diamonds) sources. 

At larger quark mass the pseudoscalar-scalar difference degrades as a pseudoscalar source, compared to pseudoscalar 
or axial currents. The window of timeslices in which the fitted mass is stable shrinks. Presumably what is happening 
is that the heavier pseudoscalar feels the effect of the scalar meson (with which it would become degenerate in the 
heavy quark limit). This effect has been noted by Ref. |27| . 

Note (parenthetically) that the data strongly hint at the presence of an increase in mp s /m q at small m q which is 
consistent with a quenched chiral logarithm. The data set is not optimized to precisely pin down these effects (the 
smallest quark mass is still rather large). In an attempt to look for these terms, I have fit m^/m q to two functional 
forms: The first one is [igll2fjj: 

(am PS ) 2 /(am q ) = A[l - 5(ln(Am q /A) + 1)] + B(am q ) (15) 

It contains three free parameters, A, B, and S, as well as the scale for the chiral logarithm, A, which I do not let float 
in the fit but instead pin to any of several fiducial values between 0.8 and 1.2 GeV. 
The second functional form is due to Sharpe [Tsj : 

(am PS ) 2 /(am q ) = C(am q ) (s/(1+S) + D(am q ) (16) 

I elect to perform fits to a mix of pseudoscalar-scalar difference correlators at low quark mass and axial correlators at 
high mass (rejecting the pure pseudoscalars on theoretical grounds). I saw little shift in the fitted S as I varied A, and 
the fits give 5 = 0.17(5) at [3 = 5.9, 0.20(7) at = 6.1. The power law fit gives S = 0.18(5) at j3 = 5.9 and 0.21(9) at 



9 



(3 = 6.1. (All chiral extrapolations from the (3 = 6.1 data set are noisier than their j3 = 5.9 counterparts because the 
former set does not extend to as low a value of the quark mass.) 

Recent determinations of S from a variety of actions show quite wide spread (for a recent review, see '28\). I am 
very much aware that my results have uninterestingly large error bars compared to the recent high-statistics, low 
pseudoscalar mass results of Ref. |2^|. My purpose in quoting them here is to illustrate that I see behavior in the 
data which is consistent with the expectations of quenched QCD, and inconsistent with the expectations of full QCD. 
The choices for the quark mass relevant to the kaon will occur in the rise near the minimum of the m 2 PS /m q vs m q 
Plot. 

These results are quite a bit greater than ones reported by Heller and me [3(|, using an overlap action with the 
same fermion kernel, but with gauge connections which are much more smoothed than the ones used in this study. 
Presumably the smaller value for 5 seen there came about because the additional smoothing suppressed the coupling 
of the fermions to smaller topological objects. 

The vector meson mass is shown as a function of the quark mass (both in lattice units) in Fig. [5J. It was measured 
with the correlator of two vector currents. It is expected to have only an "interference term" like the axial vector 
correlator. As the quality of the signal degrades at low quark mass, and as I only use vector masses at larger quark 
masses, or use all masses to make linear extrapolations to the chiral limit, I will not investigate this data set further. 
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FIG. 5: amy vs am q , (a) /3 = 5.9, (b) (3 = 6.1. 

The pseudoscalar decay constant, defined through fpg = mp s /(2m q )(0\'tlij 5 ip\PS} , is shown in Fig. EJ It is 
extracted it from a correlated fit to a two-point function with Gaussian source and sink and one with a Gaussian 
source and local pseudoscalar current sink. I have used pseudoscalar and axial sources, as well as a fit where the 
two correlators are the pseudoscalar-scalar difference. There is some tendency for the latter correlators to produce 
a smaller decay constant at smaller quark mass and a larger one at larger quark mass. The deviation of the decay 
constant measured using pseudoscalar sources and sinks from that measured in other channels at lower quark mass 
may reflect the effects of zero modes in that channel (it is consistent with previous work by[3l]]). The deviation of the 
quantity measured in the pseudoscalar-scalar difference channel at larger quark mass is almost certainly due to the 
degradation in the signal from the difference correlator there: the fit simultaneously produces the pseudoscalar mass 
along with the decay constant, and this mass also drifts away from that of the other channels. 

When (afps)is needed below, I will use the pseudoscalar-scalar difference at low quark mass and the axial source 
correlators at higher mass, and vary the crossover mass to make sure that the results are not sensitive to it. 

Finally, there are the operators contributing to Bk itself. I performed correlated fits of three and two point functions 
(Eqs. 1131 and l*HE|l to extract the matrix element of 0+. I considered two cases: pseudoscalar sources for the three 
point function, and the pseudoscalar - scalar difference for the two point function, and fits to correlators where all 
interpolating fields are axial currents. Differences in the fits under variation in the ranges are comparable to the error 
bars shown. I also investigated fits to the axial vector current operator, using source interpolating fields which were 
either pseudoscalar or axial vector. An example of such a search is shown in Fig. Q. Only at the smallest mass, 
not at the kaon, does the operator appear to be affected. Thus the extraction of the B parameter in the chiral limit 
might be compromised by finite volume effects, but not Bk itself. 
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FIG. 6: (a) (afps) vs m q , (a) f3 = 5.9, (b) (3 = 6.1. The labels crosses, squares, and diamonds label decay constants extracted 
from pseudoscalar correlators, differences of pseudoscalar and scalar correlators, and correlators with axial sources, (c), (d): 
The small mass region of panels (a), (b) magnified. 



IV. DATA ANALYSIS FOR B K 
A. Fits to Lattice Bk 



Bk can be extracted from covariant fits or ratios of three propagators. A typical correlated fit of the figure-eight 
graph and two axial vector sink correlators, and a ratio plot are is shown in Fig. [SJ Correlated fits are stable against 
variations in fit ranges or cuts on the correlation matrix and have reasonable confidence levels. Uncorrelated jackknife 
"ratio fits" are consistent with the correlated fits, have small uncertainties and are also quite stable over a wide range 
of time slices. Examples of these extractions are shown in Figs. El EH arL d ^2 My results for the B-parameter in 
lattice regularization as a function of the bare quark mass are shown in Tables III and ITT1 
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FIG. 7: Three parameter correlated fits extracting the matrix element (0+) from a correlator of the operator and a two point 
function, (a) /3 = 5.9, (b), f3 = 6.1. Meson interpolating fields are squares for pseudoscalar-scalar difference source and octagons 
for the axial vector source for both correlation functions. The "denominator" of Bk can be checked by extracting the axial 
vector current matrix element mpsfp S from a correlated fit of two point functions; this is shown in (c) and (d) for j3 = 5.9 and 
6.1. Here crosses show the pseudoscalar-axial vector correlator and octagons from the difference of (7075 — 7075) and (70 — 70) 
correlators, which has no zero modes. 



B. Matching Factors 



The naive dimensional regularization (NDR) Bk at a scale \i is related to the lattice-regulated number (from lattice 
scale a) by 

- RB K att) = ZTOT }^ a) B K att) (17) 

where the conversion factor for the operator + , Ztot, can be decomposed into 

Z T ot{h, a ) = Z NDR (n, m)Z PT (ma, qa). (18) 

Za is the axial current matching factor. Z NDR (fi,a) is the usual two-loop running formula for the NDR operator, 
carrying it from scale m to scale fi. ZpT(ma, qa) is the lattice-to-NDR matching factor connecting lattice scale a with 
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FIG. 9: Bk extracted from (a) correlated fits to three propagators and (b) jackknifed ratios of propagators, with pseudoscalar 
sources, (3 = 5.9. Fits are over the range t = 6 — 10, 24-28. In (a), squares show three-parameter fits and crosses show 
four-parameter fits, while in (b) the squares are from fits of the ratio to a constant, while the crosses include the lowest p / 
contribution. 



continuum scale m, 

Z 'prima, qa) = 1 + %. J [— 7o In ma + b] (19) 

lD7T z 

7o is the anomalous dimension of the operator; b is the lattice-action-dependent matching term (b = —3.98 for the 
action used here 0])- One can imagine several different choices for the product of Z's. 

• ma = 1, q = q* as defined by the Lepage-Mackenzie matching condition ^(|. I will call this "type 1." 

• Gupta, Bhattacharya, and SharpeQ suggest two versions of "horizontal matching," where m = qi, q* = q{. and 
qia = 1 (type 2) or ir (type 3). 
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FIG. 11: Bk extracted from jackknife averages of ratios of correlators, with axial vector sources, (a) f3 = 5.9, (b) ft — 6.1. 
Squares show constant-ratio fits and crosses show fits keeping the p^0 term. 



The RGI (renormalization group invariant) Bk is computed from the NDR B-parameter in the usual way, and I will 
write the constant which renormalizes the lattice number as Bk = RB K att ^ . 

I determine the strong coupling constant (in the so-called ay prescription) at a scale q — 3.41/a from the plaquette, 
convert it to MS prescription, and run it to the desired scale using the two-loop evolution equation. (In all cases I 
determine Z A perturbatively, 0.988 at p = 5.9, 0.990 at 6.1.) At /3 = 5.9 the "type 1" Z-factor is R = 0.97 and the 
type 2 choices are 1.00 and 0.95. The conversion factors to the RGI B are 1.35, 1.40, 1.31 for the three choices. The 
/3 = 6.1 numbers are almost equal: for B K VDR) (n = 2 GeV), 0.98, 1.01,0.97 and for B, 1.38, 1.42, 1.36. It appears 
that a reasonable choice and uncertainty for R is 0.97(3) at P = 5.9 and 0.98(3) at 6.1. (Note that in an extrapolation 
to the continuum, the same choice for matching convention must be used at all lattice spacings.) 
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C. "Quenched phenomenology" and the strange quark mass 

Bk comes from a meson made of two degenerate quarks with a common mass which is the average of the nonstrange 
and strange quark mass. Unfortunately, if quenched spectroscopy is different from the spectroscopy of full QCD, there 
is no unique choice for this quark mass, even in the absence of discretizaton-induced scale violations. What follows is 
a discussion which lies outside the realm of lattice calculations, as I try to motivate some choices for lattice spacings 
and quark masses over others. Unfortunately, Figs I9I1UI show a rather large dependence of the B parameter on the 
quark mass at small quark mass, so the choice of mass is important phenomenologically. 

One choice would be to take the lattice spacing from the Sommer parameter (using the interpolation formula of 
Ref. 32] ) and to set the pseudoscalar mass at the kaon. Either interpolation or a quadratic fit to the pseudoscalar 
mass mp s = dm q to ttiktq = 1.256 gives am s /2 = 0.027(1) at (3 = 5.9, 0.019(1) at 6.1. (The reader might recall that 
the interpolation formula of Ref. |32j ] plus a nominal value of ro = 0.5 fm implies that a would be 0.11 or 0.08 fm at 
(3 = 5.9 or 6.1.) 

However, setting the lattice spacing from ro will push the predictions for other physical quantities away from 
experiment. For example, one can extrapolate the vector meson mass, the pseudoscalar decay constant, and their 
ratio, into zero quark mass. (Both extrapolations are linear in m q ; in contrast to full QCD, the quenched decay 
constant has a vanishing chiral logarithm.) A jackknife extrapolation of fps/m p gives 0.166(6) at f3 = 5.9 and 
0.169(8) at (3 = 6.1, quite close to the physical value of (132 MeV)/(770 MeV)=0.17. Separate extrapolations of my 
and fps to zero quark mass would combine with the Sommer parameter to produce too small values of the rho mass 
or pseudoscalar decay constants: jackknife extrapolations of the former quantity give am p = 0.51(2) at f3 = 5.9 and 
0.35(1)) at (3 = 6.1, while afps = 0.084(2) or 0.060(2). This means that a better match to hou-Bk phenomenology 
requires using either f n or m p to set the scale. 

I tried a variety of fits and dimensionless ratios to find a quark mass for the kaon. Two choices which produce 
similar (but not identical) values are (1) Take the lattice spacing from the rho mass; interpolate or fit m PS vs m q 
over a range of quark masses around the kaon to the kaon mass. (2) Simultaneously determine a and m s from the 
phi meson (assume linear variation of the vector meson mass on the quark mass and the "strange eta" ss meson: 
mi + m 2 = 2m| f ); then shift the quark mass: (m s + m ns )/m s = 2rn? K /ml e . Both fits are taken over a range of 
quark masses between l/2m s and 3/2m s . One can substitute polynomial interpolation for fitting with no change in 
the resulting quark mass. The first choice gives am q = 0.035(2) at f3 = 5.9, 0.025(2) at (3 = 6.1; the second choice, 
am q = 0.029(2) at /3 = 5.9, 0.023(2) at /3 = 6.1. The values of the /j, = 2 GeV MS strange quark masses from these 
two choices 110(10) or 92(10) MeV at (3 = 5.9 and 112(7) or 104(8) MeV at = 6.1. There is no way to determine if 
the difference in these numbers is statistical or a quenching systematic. 

A comparison of B^ dr (/j, = 2 GeV) as a function of the pseudoscalar mass in units of the kaon mass for this choice 
of quark mass ((mps /itik) 2 = rn q /((m s + mj)/2)) is shown in Fig. 1121 The two data sets appear to scale reasonably 
well. 

In practice, I perform the interpolation of the B-parameter to the kaon under the same jackknife as the determination 
of the strange quark mass. Both determinations have a statistical uncertainty, and the jackknife attempts to include 
their correlations. I tried two methods of interpolation: I fit the data of the lowest five quark masses to the chiral 
logarithm form (Eq. 1201 below) or just did a simple polynomial interpolation. Both results agree within uncertainties, 
although at (3 — 6.1 the results of the chiral fit produce unacceptably large uncertainties because the chiral parameters 
are poorly determined. I will use the interpolated values. 

To be definite, make choice (2) for the quark mass and take the Lepage-Mackenzie choice for the Z-factor. In that 
case, one will have B^ DI H^ = 2 GeV)= 0.540(25) at f3 = 5.9, 0.546(36) at /? = 6.1. They lie a bit lower than the 
staggered JLQCD resultJM or the other overlap calculation of Garron, et al0], but bracket the two domain wall 
results of Refs. ^(| an d UM- The alternative "choice (1) quark mass" would give Bk = 0.563(23), 0.552(28) for the 
two couplings, 

Recall that Ref. se ^ s a from the rho mass and chooses the bare quark mass to put tuk at "its physical mass" 
while Ref. ^(j sets the strange quark mass "from the experimental value of m^/m p ." My neglect of ro as compared 
to hadronic masses is at least consistent with the choice of groups with which I want to compare. (Of course, that does 
not make it correct.) Oddly enough, the "choice 1 mass" which gives us a larger Bk is more like the RBC convention 
(and their value for Bk is below mine), and the "choice 2 quark mass" which gives a lower mass is more like that of 
Ref. [T^j, which lies above my data. Perhaps I am just being obsessive about one standard deviation effects. 

Although they show no scaling violations, an extrapolation of the "type 2 mass" results to zero lattice spacing 
B(a) = Bq + a 2 B\ would give a continuum prediction of Bk = 0.550(67) for this choice of quark mass and matching 
factor. (The RGI B K would be 0.751(35) at 5.9, 0.769(51) at 6.1 and 0.781(94) for an extrapolated continuum 
value.) These values are the points shown in Figs. ^ and [21 The alternative "choice (1) quark mass" would give 
Bk = 0.544(53) for the continuum value, and any other of the choices for the perturbative Z factor would produce 
a continuum result whose central value is also within 0.01 of Bk = 0.55 with an extrapolation uncertainty of 0.05 
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FIG. 12: B% DR (v = 2 GeV) as a function of quark mass from the two data sets-diamonds for /3 = 5.9, squares for /3 = 6.1. 
The matching factor uses the Lepage-Mackenzie convention and the quark mass comes from a combined fit to the cf> and ss 
pseudoscalar. 

(Bx = 0.78(8)). The variation in Bx at fixed lattice spacing from varying the choice of Z is completely washed out 
by the extrapolation. 

These results are quite similar to continuum 1/N C predictions of B = 0.7(1) [33| . 



V. RESULTS IN THE CHIRAL LIMIT: 0+, Or AND O s 

A. 0+ 

The behavior of Bk as a function of pseudoscalar meson mass has been computed by Sharpe[l9|. His result 
(simplified to the degenerate-mass limit appropriate to these simulations) is 



B(m)=B[l-^^\ogm 2 



87T2/2 



PSJ 



l PS 



(20) 



I follow the path of Ref. [Tl| and perform a chiral extrapolation using the functional form of Eq. 1201 The range of 
validity in quark mass of this formula is not known a priori, and so I simply performed fits of my data to it dropping 
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one heavy quark mass at a time, and looking for stability in the extrapolated result. The fit involves a jackknifc 
to incorporate the measured pseudoscalar mass, B-parameter, and /„.. A "typical" set of fits (one of the ones in a 
jackknife ensemble) is shown in Fig. EH I find that the quality of the fit is good, and extrapolated value of B is 
insensitive to the number of masses kept, as long as I keep less than about 8 of the masses (this corresponds to a 
pseudoscalar/ vector mass ratio < 0.75) I find extrapolated "type 1" values of B = 0.27(2) at (3 — 5.9; 0.23(4) at 
f3 = 6.1 (or B = 0.36(2), 0.31(7)). This is in reasonably good agreement with the more precise RBC number of 
B = 0.267(14). 




FIG. 13: A set of chiral extrapolations over various mass ranges (shown by the extent of the solid lines) to lattice Bk, from 
jackknifed ratios of correlators at /3 = 5.9 . 

The authors of Ref. performed a chiral extrapolation of their data letting the coefficient of the chiral logarithm 
be a free parameter. Their result (0.420(16)) is much larger than mine. I believe that this is the result of using 
a different functional form to extrapolate. I believe that using the known analytic functional formula for chiral 
extrapolation is theoretically better justified. 



B. Oj and 8 

My conventions for the operators and matrix elements are the same as those of Refs. and [25|, and have been 
given in Eqs. and [7| In the degenerate quark mass limit where I work, matrix elements of these operators only 
involve figure-eight graphs. The physical parameters one wishes to associate with these operators are the zero-quark 
mass limit of the matrix elements, and so a chiral extrapolation is needed. O7 and 0% also mix with each other under 
renormalization. The mixing matrix was computed in one loop perturbation theory in Ref. |l5j| . I compute the NDR 
matrix element by extrapolating each lattice operator to the chiral limit, then mixing and running the operators, all 
inside a single elimination jackknife. 

The matrix element of these operators has a chiral extrapolation similar to that of Bk, 

< n\O t \K >= C(l + Hm% s )) + bm 2 PS . (21) 
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At the time that Ref. |Tl| appeared, the size of the chiral logarithm for these operators was not known, and the 
authors did fits which either dropped the logarithm or left it as a free parameter. Since then, Golterman and Pallentc 
[34l | have computed the coefficients for chiral logarithms in quenched and partially-quenched QCD, and found that 
£ = for O7 and Og (in the degenerate- mass limit used here). Accordingly, I just extrapolate the data to zero 
assuming a linear dependence on the quark mass. 

A peculiarity of the quenched approximation is that the B-parameters for these operators go to zero in the chiral 
limit, because m 2 p S /m q diverges. Recall that the B-parameters for these operators are defined as 

B ^ = . . {K^awj (22) 

c P (7r|<075^|O) (0\^ 5 ip\K) + c a (^|^7m75^|0) (0| ^7^75^1^) 

where the Cj's are numerical coefficients. The PCAC relation says that (O^jstplPS) = \mp s /m q fps- The divergence 
of (mpg/mq) 2 in the denominator is not compensated by any singular behavior in the numerator. 

Therefore, I focus on the matrix elements themselves. Several forms have been presented in the literature; {K\Oi\ix) 
itself, 

{K\OM) = 12^ (23) 

Jtt 

the lowest order (in chiral perturbation theory) K — > (7tit)i=2 matrix element 

M i{ K - Ml=2) = ffi^ (24) 

Jtt 

and a dimensionless coupling constant, defined as 

(™), =2 ) = (25) 

Jtt 

Note Oi = 0^"(l/a) 4 , a 4 = a^ att (l/a) 6 , and M, = M\ att (1 / a) 3 . 

Because one is interested in the matrix elements themselves, I extract them from correlated fits to a three-point 
(figure-eight) function and a two point function with the same source and sink. As before, I have measured three 
point correlators with pseudoscalar and axial vector Gaussian sources. I expect these operators might show different 
zero mode effects. In Fig. ^]I show (O7) and (Og) for two choices of source/sink (pseudoscalar and axial) for the 
three point function, and pseudoscalar sources and sinks, the pseudoscalar-scalar difference for the two point function, 
and axial source and sink for both correlator. Only at the lightest quark masses are there noticible differences. I have 
performed jackknife extrapolations of the matrix elements to the chiral limit, using the lightest five quark masses 
and taking a simple linear mass dependence. When I do that, I find that the uncertainty on the extrapolated matrix 
elements inflates so much, that no dependence on the choice of sources and sinks is seen. 

I extrapolate the alternative versions of the matrix elements by including the calculation of the pseudoscalar decay 
constant in the jackknife and extrapolation. (For example, the Mi's are found by extrapolating (Oi(m q )) / fps(m q ).) 
The matrix elements and decay constants are both gently falling functions of the quark mass, and their values 
are strongly correlated in the data set. Thus, different combinations will extrapolate differently. My experience is 
illustrated with the noisier /3 = 6.1 data set (Fig. I15|) : The Oj's vary strongly with quark mass and the error on the 
extrapolated value is large. Ratios of (Oi(m q )) to powers of fps fare better. However, the high dimensionality of 
the operator, Mj, and a, mean that my results for these operators are particularly sensitive to the choice of lattice 
spacing. Thus I will restrict myself to a presentation of the matrix elements themselves and the dimensionless <7i's, 
shown in tables ITTllIVI 

The three choices of matching factors correspond to the following mixing matrices at p = 5.9, 6.1: Type (1): 

„_/ 0.979 0.005 \ / 0.979 -0.060 \ , , 

' -0.027 1.099 ) ' I -0.019 1.018 ' 



Type (2): 



Type (3): 



0.979 0.0016 \ / 0.996 -0.040 
-0.038 1.110 J ' ^ -0.041 1.000 



(27) 



/ 0.996 -0.026 \ / 1.010 -0.059 \ , , 

^-0.042 1.021 J \-0.0U 0.934 J ^ ' 

Since the type 1 and 2 Z-factors are so similar, there is a negligible difference in the results of the matrix elements 
from using either. Using the Type 3 Z's suppresses the operators by ten per cent at P = 5.9, five per cent at p = 6.1. 
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FIG. 14: (a) (J3 = 5.9) and (b) (J3 = 6.1) (0 7 ) and (c) (fi = 5.9) and (d) (J3 = 6.1) (Og) from (diamonds) pseudoscalar sources 
and sinks in the three point function and the pseudoscalar-scalar difference for the two point function, and (cross) axial source 
and sink for both correlators. 



VI. CONCLUSIONS 



Using a lattice action with exact chiral symmetry to do a lattice calculation of a matrix element in quenched 
approximation is a two edged sword. On the one hand, exact chiral symmetry eliminates "non-continuum" operator 
mixing and makes the calculation reasonably straightforward. On the other hand, because one can do simulations at 
small values of the valence quark mass, the artifacts of the quenched approximation are clearly revealed. It is clear 
that quenched QCD is not even qualitatively the same theory as full QCD, and that using the quenched approximation 
to compute hadronic matrix elements is merely another kind of phenomenology. 

My results show no scale violations. Their continuum extrapolation is the prediction for quenched Bj^ dr \/j, = 2 
GeV) = 0.55(7) or B = 0.79(9). This result brackets the two domain wall predictions of Refs. [ljj, EH, One could 
improve the statistical uncertainty of this number by further simulations, but I believe that the greatest source of 
theoretical uncertainty is an uncontrolled systematic error arising from the use of the quenched approximation itself. 
I know no way to reliably quantify this systematic. It would be a far better use of computer resources to devote them 
to simulations done in full QCD. 
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FIG. 15: Various scaled versions of Or and 0$ from the j3 — 6.1 data set. The data shown at the origin are the extrapolated 
values (a) (K\Or\ir), scaled by 3000 (octagons); ar, from Eq. 1231 scaled by 10 6 (crosses); Mr, from Eq. 1241 scaled by 30 
(diamonds); gr, from Eq. 1251 scaled by 0.01 (squares), (b) {K\Os\tt), scaled by 100 (octagons); as, from Eq. 1231 scaled by 
10 6 (crosses); Ms, from Eq. 1241 scaled by 30 (diamonds); gs, from Eq. 1251 scaled by 0.005 (squares). 
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APPENDIX A: CHECKING PERTURBATION THEORY 



It is reasonably straightforward to make a nonperturbative determination of the axial current matching factor Z 4 
and the lattice-to- MS matching factor for the quark mass, Zm, and check the perturbative calculation of Ref. [lFjJ. 

I determine Za as follows: The matrix element of the pseudoscalar current gives (0\4>j5ip\PS) = fp S mp s /(2m q ), 
with no lattice-to-continuum renormalization factor for fpg] fp S = fps- The zeroth component of the axial current 
has as its matrix element (0\4!jo75i J \PS) = fps m PS, with fp S = ZaJps- Thus Za = fps/fps- ^ compute the ratio 
by fitting a correlator with a space-summed pseudoscalar sink and a correlator with an axial vector sink, taking the 
ratio, extrapolating to the chiral limit, and jackknife-averaging. The result, which is quite insensitive to the range 
of quark masses or the choice of timeslices kept in the fit, is 1.00(2) at (3 — 5.9, 1.02(2) at 6.1. The perturbative 
calculation of Ref. [HI predicts Z A = 0.973 or 0.993 at /3 = 5.9, 0.980 or 0.988 at 6.1 (depending on the choice of 
lowest-order or higher-order q*; the latter numbers are "preferred" by the ideology of choice of momentum scale for 
the strong coupling constant of Ref. 0D- %A f° r P = 5-9 and 6.1 is illustrated in Fig. 

The other renormalization constant which has an "easy" comparison to a nonperturbative calculation is the matching 
coefficient for the MS quark mass, as described by Ref. 35] . The method involves determining the value of bare quark 
mass at which the pseudoscalar mass takes on a certain value and combining that bare mass with an appropriately- 
rescaled renormalization factor computed using Wilson fermions j3|| , and called U m by the authors of Ref. [37] . The 
Z— factor which converts the lattice quark mass to the RGI (renormalization group invariant) quark mass is 

Zm{9q) = U m \(r mps) 2 =x rs f (Al) 

The RGI condensate can then be converted to the MS regulated condensate using a table of (multi-loop) conversion 
coefficients from Ref. [3^: Zjjg(ji) = Zmz(^x) where z = 0.72076 for /1 = 2 GeV. 
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FIG. 16: Inferred Za from ratio of axial vector to pseudoscalar matrix elements, (a) = 5.9 (b) = 6.1. The black square is 
the perturbative prediction. The open square is the jackknife extrapolation to the massless limit. 



The authors of Ref. [26J quote U at x = rgTOp S =5, 3, 1.5376. The larger x values occur at quark masses for which 
m 2 PS /m q is a gently falling function (compare Fig. 0]and Table IV|l . For these values, polynomial interpolation of the 
pion mass gives the values shown in the Table. However, the smallest value of x corresponds to a quark mass deep in 
the region where m 2 PS /m q is rapidly rising. To drive the pseudoscalar mass to the fiducial value of x requires a quark 
mass which is "smaller than it should be" (in the absence of a rise), and so the value of Zm is enhanced, to about 
1.21(8) for the (3 = 5.9 data set. 

Note that the authors of Ref. [3(| actually did not measure U at x — 1.5376; instead, they extrapolated their 
pseudoscalar masses down from higher values of the quark mass. The extrapolation is done assuming that vn 2 PS is 
linear in m q , that is, without chiral enhancements. Simply ignoring the lower quark mass data, fit the higher quark 
mass data (am q = 0.035 to 0.125 at (3 = 5.9) to a 2 mp s = Cam q , extrapolate to lower quark mass, and using the 
extrapolated mass in Eq. IA1I produces the results shown in the Table. All three Z— factors seem to be consistent, 
and are slightly higher than the perturbative prediction of 1.00 at /3 = 5.9 (the constant term nearly cancels the 
logarithm there) and 1.02 at (3 = 6.1. 
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TABLE I: Pseudoscalar mass and lattice-regulated Bk as a function of bare quark mass, from correlated fits, at (3 = 5.9. 



arriq 


amps 


lattice Bk 


0.015 


0.227(5) 


0.476(33) 


0.020 


0.254(5) 


0.512(25) 


0.025 


0.278(5) 


0.540(27) 


0.035 


0.320(5) 


0.581(23) 


0.050 


0.386(4) 


0.636(18) 


0.075 


0.467(4) 


0.667(36) 


0.100 


0.537(3) 


0.683(36) 


0.125 


0.603(3) 


0.724(43) 


0.150 


0.665(3) 


0.753(29) 


0.200 


0.784(2) 


0.794(37) 



TABLE II: Pseudoscalar mass and lattice-regulated Bk as a function of bare quark mass, from correlated fits, at (3 = 6.1. 



am q 


amps 


lattice Bk 


0.015 


0.192(5) 


0.512(38) 


0.020 


0.218(5) 


0.535(29) 


0.025 


0.234(5) 


0.556(47) 


0.035 


0.270(4) 


0.599(45) 


0.050 


0.319(3) 


0.657(12) 


0.070 


0.379(3) 


0.720(10) 


0.090 


0.435(2) 


0.771(17) 


0.110 


0.488(2) 


0.805(21) 


0.150 


0.587(2) 


0.793(50) 



TABLE III: (3 = 5.9 atrix elements of operators O7 and Og: first matrix elements (K\Oi\n), in units of a 4 , i = 7, 8, then 
dimensionless NDR matrix elements </;. All fits use figure-8 graphs with pseudoscalar sources and "type 1" Z- factors. 



am q 


{K\Or\n} x 10 3 


(K\O s \n) x 10 2 


97 


38 


0.015 


6.83(88) 


2.14(28) 


169(38) 530(119) 


0.020 


6.00(59) 


1.86(18) 


111(17) 


343(51) 


0.025 


5.63(82) 


1.72(37) 


87(14) 


266(64) 


0.035 


5.33(66) 


1.60(20) 


67(10) 


199(29) 


0.050 


5.79(56) 


1.68(17) 


57(6) 


164(18) 


NDR extrap 


6.17(90) 


2.17(33) 


127(19) 


433(77) 
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TABLE IV: (3 = 6.1 atrix elements of operators O-j and Og- first matrix elements {K\Oi\ir}, in units of a 4 , i = 7, 8, then 
dimensionless NDR matrix elements </;. All fits use figure-8 graphs with pseudoscalar sources and "type 1" Z- factors. 



am q 


(K\0 T \n) x 10 3 


(K\O s \tt) x 10 3 


97 


58 


0.015 


1.87(33) 


5.64(108) 


197(56) 594(176) 


0.020 


1.69(26) 


5.20(91) 


125(30) 


384(98) 


0.025 


1.64(26) 


4.88(71) 


92(17) 


274(57) 


0.035 


1.82(19) 


5.30(58) 


71(9) 


208(28) 


0.050 


2.18(15) 


6.11(44) 


60(5) 


168 15) 


NDR extrap 


1.08(37) 


4.04(120) 


111(22) 


379(82) 



TABLE V: Nonperturbative determination of Zm 



X 


U 


rorriq 


Zm 


Z^(n = 2 GeV) 


13 = 5.9 










5 


0.580(12) 


0.380(5) 


1.52(4) 


1.10(3) 


3 


0.349(9) 


0.216(5) 


1.62(6) 


1.16(4) 


1.5376 


0.181(6) 


0.121(1) 


1.50(5) 


1.08(4) 


(3 = 6.1 










5 


0.580(12) 


0.385(6) 


1.51(4) 


1.09(3) 


3 


0.349(9) 


0.228(6) 


1.53(5) 


1.10(4) 


1.5376 


0.181(6) 


0.116(3) 


1.56(6) 


1.12(4) 



